Width-amplitude relation of Bernstein-Greene-Kruskal solitary waves 



o 
o 

(N 
•*-> 
o 

O 

m 



Li- Jen Chen, 1 David J. Thouless, 2 and Jian-Ming Tang 1 

department of Physics and Astronomy, University of Iowa, Iowa City, IA 52242-1479 
2 Department of Physics, University of Washington, Seattle, WA 98195-1560 
(Dated: February 2, 2008) 

Inequality width-amplitude relations for three-dimensional Bernstein-Greene-Kruskal solitary 
waves are derived for magnetized plasmas. Criteria for neglecting effects of nonzero cyclotron radius 
are obtained. We emphasize that the form of the solitary potential is not tightly constrained, and 
the amplitude and widths of the potential are constrained by inequalities. The existence of a con- 
tinuous range of allowed sizes and shapes for these waves makes them easily accessible. We propose 
that these solitary waves can be spontaneously generated in turbulence or thermal fluctuations. We 
expect that the high excitation probability of these waves should alter the bulk properties of the 
plasma medium such as electrical resistivity and thermal conductivity. 
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Coherent structures with nonuniform charge densities 
are ubiquitous in plasma systems. Laboratory experi- 
ments have shown that such structures can be generated 
by applyingvoltage pulses Q, 0] , voltage jumps , in- 
tense laser or plasma beam injections 5]. Increasing 
numbers of space-borne observations have revealed fre- 
quent appearance of electrostatic solitary structures in 
space plasmas (0, and references therein) including 
regions where magnetic reconnection occurs 0, Q . Soli- 
tary waves can efficiently transport energy, momentum 
and charge, and are one of the building blocks in a de- 
terministic description of turbulence |lOt . The study 
of their allowed parameter space is crucial in establishing 
their relevance to real systems. Most solitary waves, such 
as those for shallow water (Korteweg-de Vries solitons) 
and those that describe crystal dislocations (sine-Gordon 
solitons) , have a strict one-to-one map ping; between their 
widths, amplitudes and velocities |12j . However, in colli- 
sionless plasmas, the widths and amplitudes of the elec- 
trostatic solitary waves that exhibit vortex structures in 
phase space |l3( are not tightly constrained [lj, [l5j . Al- 
though the width-amplitude relation of these so-called 
Bernstein-Greene-Kruskal (BGK) electron solitary waves 
in three dimensions (3D) have been studied by a few au- 
thors [bH Hjl H(i ] , the results were limited either by im- 
complete analysis of the width-amplitude relation [16| or 
by incomplete solutions |14l Il5| , and the broad impact 
has not been recognized. In this letter, we derive in- 
equality width-amplitude relation for 3D BGK solitary 
waves. The description is valid for both electron and ion 
modes. The inequality width-amplitude relation dictates 
a continuous range of allowed sizes and amplitudes for 
these waves, and thus makes them easily accessible. We 
propose that the continua of allowed existence range en- 
able BGK solitary waves to be spontaneously generated 
in thermal fluctuations as well as turbulence, and is re- 
sponsible for their ubiquitous presence in widely different 
classes of collisionless plasmas. 

To construct exact nonlinear solutions that are local- 



ized in 3D, we use the BGK approach that was for- 
mulated for ID nonlinear Vlasov-Poisson equations [l7j . 
but extend the Poisson equation to 3D. We construct 
azimuthally symmetric solutions in the limit of infinite 
magnetic field, and then tune down the magnetic field to 
obtain the criteria for neglecting effects of nonzero cy- 
clotron radius. One key step in the BGK approach is 
to separate particles that are trapped in the potential 
and those that are passing. We prescribe the potential 
form and the passing particle distribution, solve for the 
trapped particle distribution, and derive the physical pa- 
rameter range. This approach is much easier than to 
prescribe the passing and trapped particle distributions 
to solve for the potential [131 Il8| . and so allows us to 
explore the solution space much more fully. 

We consider two species of charge carriers (electrons 
and one type of ions), each with charge q Sl mass m s , and 
thermal energy T s . The background magnetic field B is 
along the z direction. In the strong field limit, particles 
move only along B with velocity v, and the distribution 
functions f s satisfy the following Vlasov equations, 



df s (r,v) q s <9$(r) df s (r, v) 



dz 



dz 



dv 
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where $(r) is the electrostatic potential. It can be easily 
shown by chain rules that any f(r±,w) is a solution to 
Eq. J3> since its dependence on z and v is only through 
the particle energy w = m s v 2 /2 + q s &(r). Such distribu- 
tion functions and the potential are further constrained 
by the Poisson equation, 



V 2 $(r) = £ / 



dw 4im, g ,/,(r ± ,u>)_ (2) 
i-Az s *(r) y/2m s [w - q s ®(r)] 

where the velocity space integral has been converted to 
energy space integral, and f s (r±, w) is normalized so that 
n s is the particle density in the unperturbed region where 
charge neutrality gives J2 s ns< l s = 0- Species 1 is defined 
to be the one which involves trapping (min [gi$(r)] < 0), 
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and species 2 does not. The distribution function fi is 
further divided into passing and trapped components, f p 
and f tr - The second term in Eq. JTJ is nonlinear as <I> is a 
functional of the particle distributions and vice versa. In 
physical terms, the system is nonlinear because plasma 
particles collectively determine the mean-field potential, 
and the potential in turn determines how particles dis- 
tribute themselves. It is the presence of this nonlinear 
term that admits solitary wave solutions which exhibit 
localized structures in potentials and distribution func- 
tions. 

Eq. can be thought as a set of ID Vlasov equations 
in the z direction for given rj_. These parallel Vlasov 
equations are coupled by the perpendicular profile of the 
potential $ through Eq. If $ is known, Eq. © re- 
duces to a set of ID integral equations parameterized by 
r_i_. For given fa and the passing distribution f p , the 
trapped distribution f tr can be found by solving these 
integral equations. The important requirement for the 
solutions to be physical is that the trapped distribution 
ftr so determined should be nonnegative. This leads to 
a self-consistent constraint on the form of the potential 



specified at the beginning. It turns out that neither the 
potential forms nor the passing distributions are tightly 
constrained. One can prescribe different localized poten- 
tial functions or different passing particle distributions 
(as long as the distribution functions satisfy the Vlasov 
equation). As an illustrating example, the solitary po- 
tential is chosen to be an azimuthally symmetric double 
Gaussian, 

&(r,z) = gi> exp (-z 2 /25 2 z - r 2 /25 2 r ) , (3) 

where g = — signal ) in order for $ to trap species 1 
particles, ip is the potential amplitude and is positive, 
r = | r_i_ | , 5 Z and 5 r are the parallel and perpendicular 
widths. The distributions f p and ji are of the Boltzmann 
type, 

f p (w) = V2m 1 /7rT 1 exp(- W /T 1 ) , (4) 
h(w) = v /2m 2 /7rT , 2exp(-w/r2) . (5) 

By carrying out the integrals of f p and /a in Eq. |J2J, we 
obtain the trapped particle density, 



n tr ($) = $ 
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where t = q2Ti/(qiT 2 ). To simplify the expression, we have set the length unit to be = \jT\j 'Airniq 2 , energy 
unit T\ and charge unit qi ($ becomes strictly negative in units of T\/qi). We have also rewritten the expression in 
terms of the potential, a crucial step for analytically solving the integral equation. We refer readers to Appendix B 
in Ref. fl9j for the details of solving this Volterra type integral equation, and simply write down the result here, 



ftrjr, w) 
y/2m 1 



r^2 



<r i 1 1 
S 2 [ft ~ 5* 
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[1 - erf (s 



= v / ierfi(v /z tw) , (7) 



where w < 0, and erfi(z) = erf(iz)/i is the complex error 
function which is a real function of its argument. 

At this juncture, it is instructive to consider the small 
amplitude behavior of our solution and compare with pre- 
vious results. Expand the RHS of Eq. ijfjj for small $. 
The leading term comes from the passing particle density 
and is of order v^, and the next order is linear in $. This 
means that for small potential amplitudes, the nonlinear- 
ity actually dominates over the linear response. In this 
case, even when Ti ~ T e for electron-proton plasma, the 
trapped distribution f tr can be solved keeping only the 
leading term, and it yields a constant •v/2mi/7r, which 
demonstrates that the plasma can sustain the solitary 
structure by providing a uniform trapped particle distri- 
bution. This result, in the context of ion holes, sharply 
contrasts the previous result 01 which predicts that ion 
holes do not exist when T e /Ti < 3.5. 



In recapitulation, we have made one step forward to 
obtain solutions in Eq. for arbitrary t taking into ac- 
count density perturbations of both species. We also ob- 
tained new understanding: nonlincarity dominates over 
linear response in the small amplitude limit - an unex- 
pected result as one normally expects nonlincarity to be 
unimportant for small amplitudes. Hereafter, we shall 
use electron holes as the example to discuss the solu- 
tion behavior and demonstrate how an inequality width- 
amplitude relation is obtained in the t — > limit, the 
limit which all previous analytical solutions for elec- 
tron holes were based on. This limit is also applicable 
to ion holes generated in laboratory experiments where 
Ti/T e - 0.03 to 0.1 . 

We illustrate below how the tuning of the parameters 
(if), S r , 5 Z ) affects the trapped particle distribution. Fig.Q] 
plots ftr (thick lines) and f p (thin lines) at r = and 
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FIG. 1: Velocity distributions of particles at the center 
(r = 0, z = 0) of BGK solitary structures for different po- 
tential amplitudes and sizes demonstrating how the tuning 
of these parameters affects particle distributions. Parameters 
are: (a) (if},8 r ,8 z ) — (1.45,5,3) (solid lines) and (ip,8 r ,8 z ) — 
(0.5,5,3) (dashed lines), (b) (ip,8 r ,5 z ) = (1.45,5,3) (solid 
lines) and (ijj,8 r ,S z ) = (1.45,5,10) (dashed lines). In both 
cases the thick lines represent trapped particle distributions, 
and the thin lines are passing particle distributions. 




FIG. 2: The inequality relation between the parallel size (8 Z ), 
the perpendicular size (<5 r ) and the potential amplitude (tp) 
showing that there is a continuum of allowed heights and 
widths for BGK solitary waves. The allowed region, the 
shaded surface and above, is marked by O and the forbid- 
den by x. The curves on the 8 Z — plane are projections of 
constant S z contours on the surface. 



z = as a function of velocity v. The solid curves in both 
(a) and (b) correspond to a BGK solitary wave that has 
zero phase space density at its phase space center (r = 
0, z = 0, v = 0). When the size of the structure is fixed, 
decreasing the amplitude raises the center phase space 
density as shown by the dashed curve in Fig. da). On 
the other hand, increasing the amplitude would lower the 
center phase space density from zero to a negative value 
(not shown), and hence result in unphysical solutions. 
When the amplitude is fixed, increasing the parallel size 
raises the center phase space density as shown in Fig.^b) 
by the dashed curve. Varying S r results in a similar effect. 

We now proceed to derive the width-amplitude in- 
equality relation. For S r < S z , the first term in Eq. J7J) is 
positive, hence the global minimum is ftrir — 0, w = —ijj) 
and its being nonnegative determines the correspond- 
ing width-amplitude relation. For 5 r > 8 Z , the global 
minimum of ft r occurs at maximum allowed r. For a 
given w < 0, the maximum r at which a trapped par- 
ticle with energy w can exist is the r max that satisfies 
— w = Q(r max ,0). Putting Eq. © into this condition, 
we obtain r^ ax = — 2S 2 . ln(— w/tp). Since the condition 
ftr(r max ,w) > guarantees f tr (r < r max ,w) > for 
S r > 5 Z , we replace r 2 in f tr (r,w) by r^ax- The global 
minimum of / tr (r max ,w) is / tr (r maX! w = hence 
the condition /t r (r max ,iD = — ip) > yields the width- 
amplitude relation for 6 r > S z . 

The conditions ft r (f = 0, w = —ip) > when 5 r < S z 
and /trf^max,!" = — VO > when S r > S z yield exactly 
the same expression. Upon re- arrangement, the resulting 
inequality is written as 

5 > / 2(41n2-l) 
z ~ y v£e*[l - erf(y^)]/v? - 4/5 r 2 " 1 ' 

Fig. [21 plots this inequality. Parameters lying on or above 
the shaded surface are allowed (O), and those under 



the surface are forbidden (x). Curves on the 5 Z = 
plane are projections of the constant S z contours to help 
visualization of the trend of the shaded surface. The 
shaded surface curves up toward infinity because the 
denominator in inequality © has to be positive, and 
that yields another relation between S r and ip (given 
by 5f > 4v / V'/(v^ r e' / '[l ~ er f VW\) an d is plotted as the 
asymptotic curve on 6 Z = plane). These inequalities 
occur because we are free to place the global minimum 
of ftr in a continuous range by correspondingly adjusting 
the amplitude and widths. A point on the shaded surface 
represents a parameter set that yields zero phase space 
density at w — —ip, that is (r = 0, z = 0, v = 0), the cen- 
ter of the solitary phase space structure. One example 
of the empty-centered distribution has been provided in 
Fig. n Lowering the amplitude or increasing the widths 
shifts a point on the surface to the region above, and 
Figs, ^a) and (b) illustrate respectively the effects on 
the distribution functions. 

We note that in the limit of 5 r — * oo, inequality JSj) re- 
duces to the width-amplitude relation for ID BGK soli- 
tary waves. This limit of inequality JSJl gives an upper 
bound for 6 Z that is valid for all finite 8 r . The result- 
ing ID inequality relation provides us a ground to un- 
derstand the discrepancy of whether the width should 
increase jlOl or decrease 0] with the amplitude. Both 
results are contained in the inequality relation with that 
by Ref. |20l | corresponding to the lower bounding curve 
since only empty-centered distributions were studied, and 
that by Ref. |l8| contained in the region above the curve 
as the distributions only take finite values in the center 
of the phase space structure. 

To establish the validity of the above results in finite 
magnetic field, we need to know whether the effect of 
nonzero cyclotron radius would result in decoherence of 
the solitary structure, that is, how the distance of the in- 
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stantaneous particle guiding center to the symmetry axis 
would vary. Therefore, it is best to look at the parti- 
cle motion projected onto the 2D x — y plane perpen- 
dicular to the magnetic field. The motion of a charged 
particle inside the solitary structure is influenced by the 
uniform B and the 3D inhomogeneous E which upon pro- 
jection becomes time-dependent E2.d- It can be shown 
that when the time variation scale of E213 is much smaller 
than the cyclotron frequency u> c , and the spatial varia- 
tion scale (L) of E2D is much larger than the cyclotron 
radius r c , the instantaneous guiding center would spiral 
around the infinite-field guiding center, and the solitary 
structure can be maintained (details of calculations will 
be presented elsewhere [ilj). Since the time variation of 
E 2 £> is characterized by the frequency (ivb) of the parallel 
trapped-particle bouncing in the potential, the condition 
can be written as 

u b /u) c < 1 — ► y/m e tjj/e/(BS z ) < 1 (9) 
r c /L < 1 — » y/2m e ip/e/(BS r ) < 1 , (10) 

where we have expressed on the right hand side of the 
arrows the condition in terms of familiar variables. For 
relevant numerical investigation of nonzero cyclotron ra- 
dius effects, the readers are referred to Ref. fl6j. 

We note that the size and the amplitude of BGK soli- 
tary waves do not have a lower cut-off within our the- 
ory. The underlying reason is that the screening of the 
charged core is accomplished by trapped particles which 
are part of the solitary structure itself. Debye screen- 
ing is not involved in these self-consistent, self-sustained 
nonlinear objects. Their size can be well below the Debye 
radius as far as there are enough particles in the solitary 
wave to ensure the validity of the mean-field approach. 
Taking a Debye radius (Ad) 100 m and a plasma density 
5 cm -3 (typical of the low altitude auroral ionosphere), a 
width of 0.01 \d for the solitary potential allows 5 x 10 6 
particles in the structure. Indeed, sub-Debye scale soli- 
tary waves have been observed |(| . 

It is the primary purpose of this paper to address the 
physical origin and significance of the inequality width- 
amplitude relation plotted in Fig. [21 The freedom to 
continuously adjust the global minimum of the trapped 
particle distribution is due to the collisionlcss nature of 
the plasma. Collisionlessness preserves the identity of 
trapped and passing particles as the energy of a particle 
is conserved. Collisions destroy energy conservation of 
individual particles, and consequently, do not allow the 
existence of trapped particle state nor the adjustment 
on the occupation number of the state. Therefore, the 
kinetic solitary waves have a continuum of allowed po- 
tential heights and widths, in great distinction to fluid 
solitons. Moreover, even with the same ambient plasma 
distribution, different functional forms for the solitary 
potential is allowed. The impact of this multitude of 
continua of allowed potentials is that these BGK states 



can be easily excited. In a system with certain fluctua- 
tion level and different fluctuation lengths, BGK states 
may be accessed easily since for a fixed amplitude there 
is a wide range of allowed widths. We therefore propose 
that in turbulent systems, BGK solitary waves can be 
spontaneously generated in the absence of two-stream or 
current-driven instabilities. The spontaneous generation 
of these coherent structures and their subsequent mu- 
tual interaction may dominate the transport properties of 
the turbulence. Because these solitary waves have many 
degrees of freedom, energy, momentum and charge are 
readily transferred between them, so that they can make 
important contributions to bulk properties of the plasma 
such as thermal transport and electrical resistivity. At 
the vicinity of the waves, passing particles are acceler- 
ated and, together with trapped particles, form coun- 
terstreaming beams (Fig. 1). Hence the velocity spread 
increases significantly and results in much higher average 
velocity spread (heating) of the plasma. Moreover, parti- 
cle trapping will prevent particles from free acceleration 
by the applied electric field and regulate the electric cur- 
rent. The high excitation probability of BGK waves can 
thus lead to finite resistivity that is required for melting 
the frozen-in magnetic flux and facilitate reconnection to 
occur in collisionless plasmas. 

In summary, we have derived inequality width- 
amplitude relation for 3D BGK solitary waves, and estab- 
lished their relevance in finite magnetic fields. The con- 
tinuum of allowed potential sizes and shapes of the waves 
is due to their kinetic nature, and leads to their ubiqui- 
tous presence. We envision that BGK solitary waves can 
be spontaneously generated in thermal fluctuations as 
well as in turbulence. 
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